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ABSTRACT 

Objectives: The objective of this paper is to analyze the unsteady hydrodynamic 
squeeze film lubrication between rough parallel plates with changed boundary 
conditions. 

Methods/Statistical analysis: The stochastic model of Christensen and Tonder has 
been deployed here to evaluate the effect of surface roughness. Also the effect of 
roughness parameters on different moments is numerically modelled. The associated 
stochastically averaged Reynolds equation is solved to obtain the pressure 
distribution. The results obtained here are presented in graphical forms. 

Findings: It has been conclusively establish that in general the transverse surface 
roughness effects the bearing performance adversely. Out of the three roughness 
parameters the standard deviation influences the bearing system significantly. Further 
a better situation is registered with reference to Xj axis. 

Application/Improvements: This theory can be developed to manufacture a 
suitable design for hip-joint. 
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1. INTRODUCTION 

Bio -tribology, the study of lubrication, wear and friction within the body has become a topic 
of encountered debilitating disease and trauma that destroy function of the joints. The hip- 
joint is a complex joint system which carries much of the weight of the body. It has 
sophisticated interaction between cartilage, bone, synovial fluid and other connected tissues. 
Due to degeneration of the joints, the use of artificial hip-joints has become necessary. With 
an increasing number of the ageing population requiring total hip replacements at a younger 
age, it is necessary to try and improve the design and longevity of the implant. A highly 
successful surgical procedure has been to replace the joint with an artificial equivalence which 
alleviates dysfunction and pain. 

The tribological performance of an artificial hip-joint has become a strong topic of 
discussion. The principle causes for failure are adverse short and long term reactions to wear 
debris and high frictional torque in the case of poor lubrication that may cause loosening of 
the implant. Thus, using experimental and theoretical approaches models have been 
developed to evaluate lubricate on performance under standardized conditions. 

In 1 studied lubrication mechanism of hip joint replacement prostheses. A laboratory study 
has been made of the lubrication and friction characteristics of a plastic on metaland a metal 
on metal total hip joint replacement prostheses. The results established the effect of lubricant, 
speed, and load on the performance of the joints. Under a 560 lb load the coefficient of 
friction of the plastic to metal joint was lower than that of the metal to metal at speeds up to 
the equivalent of fast walking. Bovine serum and synovial fluid, as well as human serum 
albumin were found to be good lubricants of both types of prostheses. The frictional force 
produced by the metal on metal prosthesis increased linearly with load, in both the dry and 
lubricated states. The results proceed that at low physiological loads the effort required to 
articulate the prostheses were comparable while at higher loads the friction force of the metal 
to plastic was significantly lower. 

In 2 analyzed the effects of proteins on the friction and lubrication of artificial joints. The 
tribological testing of artificial hip and knee joints in the laboratory has been ongoing for 
several decades. This work has been carried out in an attempt to simulate the loading and 
motion conditions applied in vivo and, therefore, the potential for the success of the joint. 
However, several different lubricants have been used in these tests. The work documented in 
Scholes et:al paper 2 compared results obtained using different lubricants and made 
suggestions for further works in this direction. 

In 3 discussed graphitic tribological layers in metal-on-metal hip replacements. Arthritis is a 
leading cause of disability, and when nonoperative methods have failed, a prosthetic implant 
is a cost-effective and clinically successful treatment. Metal-on-metal replacements are an 
attractive implant technology, a lower-wear alternative to metal-on-polyethylene devices. This 
paper dealt with the evidence for graphitic material in the tribological layer in metal-on-metal 
hip replacements retrieved from patients. As graphite is a solid lubricant, its presence helped 
to explain why these components exhibited low wear and suggested methods of improving 
their performance; simultaneously, this raised the issue of the physiological effects of 
graphitic wear debris. 

In 4 deals with the observation of model synovial fluid lubricating mechanisms. This paper 
examined the fundamental mechanisms of synovial fluid lubrication in artificial joints. The 
results indicated that two types of film were formed; a boundary layer of adsorbed protein 
molecules, which were augmented by a high-viscosity fluid film generated by hydrodynamic 
effects. The high-viscosity film was due to inlet aggregation of protein molecules forming a 
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gel which was entrained into the contact preferentially at low speeds. As the speed increased 
this gel appeared to shear thin, giving much lower lubricant film thickness. 

In 5 discussed the matter of synovial fluid lubrication: implications for metal-on-metal hip 
tribology. This which they reviewed the accepted lubrication models for artificial hips and 
presented a new concept to explain film formation with synovial fluid. The implications of 
this new mechanism for the tribological performance of new implant designed and the effect 
of patient synovial fluid properties were discussed. 

In 6 studied synovial fluid lubrication: the effect of protein interactions on adsorbed and 
lubricating Films. In this work examined how the protein content of six model synovial fluids 
affected film formation under static and rolling conditions and if the changes in properties 
could be correlated. 

In 7 viscous oscillatory fully developed flow of stratified fluid in a long vertical narrow 
rectangular channel with one side being pervious and the other side being impervious is 
considered. The outcome of this investigation reveals that the stratification effects reverse the 
flow near the center of the channel while magneto hydrodynamic effects make the flow linear 
closer to the plates. 

In 8 non dimensional equation of resistance coefficient with Reynolds number of porous 
medium flow has been studied recently. Experimental results were compared with Darcy’s 
equation and the validity of this equation was established. The derived equations were 
suggested to be applied in porous medium flow such that the velocity of flow could be 
determined from which discharge through porous medium can be estimated. 

In 9 investigated the squeezing hydrodynamics of a lubricant between two rough parallel 
plates: hip-joint representation. The stochastic model of Christensen and Tonder was deployed 
here to evaluate the effect of surface roughness. Also the effect of roughness parameters on 
different moments was numerically modelled. The results established that the standard 
deviation associated with roughness had significant impact. Further it was observed that the 
situation remained relatively better in the case of negatively skewed roughness. This effect 
advanced when variance (-ve) occurs. 

Here it has been proposed to study the hydrodynamics squeeze film lubrication between 
rough parallel plates in the context of hip joint. A new set of boundary conditions has been 
introduced. 


2. MATERIALS AND METHODS 


With usual assumptions of fluid film lubrication the velocity profiles are governed by the 
equations as presented in 10 . Besides, following boundary conditions are adopted in 10 there is 
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where h is the film thickness, the distance between surfaces. In the present situation it has 
been sought to investigate the equations 
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Conservation of mass is now taken into account by integrating equation 
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The film thickness is dependent on x x ,x 2 and t .Applying Liebnitz’s theorem of 
integration and also using boundary conditions one arrives at 

8 “'iff A , cdu , , . dh 

\u l (x x ,x 1 ,x 3 ,t)dx 3 = — -dz + u l (x x ,x 1 ,n,t) 

J J 0 dx x ' 0X| 


0X| 

a 


0 

h(xi ,x 2 ,t) 


r r c)w 

u 1 (x l ,x 1 ,x 3 ,t)dx 3 = — -dz + u 1 (x l ,x 1 ,h,t) 

"2 0 0 ^2 

But the integrals on LHS of above equations are, 

h(x l ,x 2 ,t) j J t T 7„ 2 

|m 1 (x,,X2,x 3 ,0Sx: 3 =- 


5 x n 


a/? 

ax. 


(4) 

(5) 


0 

h(x J ,^2?0 


J w 2 (Xj , x 2 , x 3 , t)dx 3 = ■ 


dp h 2 U x h cj 2 h 
dx j 12^ 2 4// 

dp h 3 U 2 h o 2 h 


dx 2 \2p 2 4^ 


( 6 ) 


In view of the stochastic model of Chirstonsen and Tonder 1 1,12,13 these equations 
transform to 
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, G(h) - h -\-3ct h + 3cc h-\-3cxh +3 cr~c(-\-(x 4-s 
where v J 

o : standard deviation, a : variance, s iskewness, h : film thickness 
Putting these values in equation (3), with the help of equation (4), (5) and resorting to 
modified boundary conditions, 
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Rewriting this relation using Kinematic boundary conditions on the surface x 3 = h one 
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where u is the surface velocity vector. 

2.1 Reynolds Equation Applied to a Hip Joint 

For modeling the Hip-joint the following assumptions are made, 

• the fluid is isoviscous (Newtonian) 

• the cup is positioned horizontally 


• walking cycle imposed is based on the Bergmann walking cycle 

Synovial fluid and blood plasma are known to have non-Newtonian fluid properties. But 
studies have concluded that at high shear rates, an isoviscous assumption is valid. The 
acetabular cup is anatomically positioned at 45° ; however, the contact mechanics model can 
be developed in a horizontal position. The Bergmann walking cycle consists of a loading 
pattern which is double-peaked. This has resulted from studies conducted by Bergmann et al. 
(1995) 14 on the influence of heel strike on the loading of the hip joint. The coordinate system 
that best describes the hip joint is the spherical coordinate system. Accordingly, expanding 
equation (9) using spherical coordinates, one comes to the conclusion that, for equation (9), 


http :// www. iaeme.com/UMET/index. asp 



editor@iaeme.com 



Hydrodynamic Squeeze Film Lubrication Between Rough Parallel Plates: Hip-Joint 
Representation with Modified Boundary Conditions 


V7D dP h 1 dP ~ 1 dP? 

VP = — R + 0 + (j) 


dR 


R dO R sin 0 d(j) 


since 

dP 


dR 


= 0 


Making use of divergence formula for spherical coordinates one concludes that 


G(h) 

12 p 


VP 


J_l_l 8_ 

R 2 sin 6 dO 


G{h ) . a dP 
-sm9 — 


12 // 


dO 


i a 

1 

f G(h) dP^] 

, sin 2 9 8(j) 

1 12// dtp) j 


For RHS of equation (9), one derives that 

V* 


'hu) 

1 

l 2 J 

2Psin# 


i nrT i ■ n dU e tt ■ tt dh 

hcosd Un +/?sin# — -+LL sin# — +h — - +£P — 

89 6 89 d(j> * d(j) 


Substituting (10) and (1 1) in (9), one arrives at 


1 


1 


R 2 sin Odd 


G(h) . .dP 
sin 9 


12 // 


dO 


sin 


2 9 d(j) 


G(h ) dP 
12 ju d(j) 


1 


2 R sin 6 


h cos 6U e +hsm6 


Sh 


DUn 

+ U B sin 0 

89 89 


, dU+ 8h 

+ h 1- U m — 

8<j) * 8c/) 


8h a 2 

+ — + ■ 


8t 4 ju 


(10) 


( 11 ) 


( 12 ) 


where the surface velocity components, U g and are given by 
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If the axes are shifted and the X x — X 3 plane lies in the horizontal position with X 2 in the 
vertical direction, the last equation can be simplified with rotation about the X 3 axis to obtain 
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The dimensionless form of the above equation turns out to equation 
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where, 

a is the standard deviation in dimensionless form, 
a is non-dimensional variance, 
a is the skewness in dimensionless form, 
p is non-dimensional pressure, 

(j) is non-dimensional rotation. 

Theassociated dimensionless boundary conditions are, 

P = 0 at ~<p=3ji 
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* Y (15) 

Integration of (13) in view of (14) leads to the expression for non- dimensional pressure 
distribution as 
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3. RESULT AND DISCUSSION 

It is easily noticed that the pressure distribution is governed by the equation (16). It depends 
upon various parameters such as a , a , s , (f> , co x and 0 . This expression makes it clear that 


a judicious choice of velocity in the X 3 direction may lead to an overall improve situation. A 
comparision of equation (16) with that of equation (15) in 9 suggest that nominally increased 
pressure is registered here due to the modification of the boundary conditions. In the light of 
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the involvement of the factor ( a + cf 2 )the graphs, which are shown in figure 1-8, are 
generated for the pressure distribution. 

It is seen that the standard deviation introduces an adverse effect by lowering pressure 
however negatively skewed roughness results in increased pressure. Same is the trend of 
pressure with respect to variance (-ve). The graphical representations indicate that the effect 

of rotation along X 3 direction remains predominant for this type of bearing system. 

4. CONCLUSION 

It is clearly observed that the modification of boundary conditions turns in a nominally 
increased pressure distribution. This study reveals that for any type of improvement, the 

standard deviation is required to be kept in minimum even if the velocity in X 3 direction is 
suitably chosen. Therefore, while designing the ball-socket bearing system, the roughness 
aspects must be given due consideration because even the negatively skewed roughness can 
increase the pressure, which may be further augmented due the occurrence of variance (-ve). 
Certainly, this study may be of some help particularly when the system has run for a long 
time. 

It was deemed appropriate to modify the boundary conditions incorporating the roughness 
parameter for deriving a better hip-joint system. Therefore, the positive effect of negatively 
skewed roughness and variance (-ve) can be channelized while designing the hip- joint 
bearing system. 



Figure 1 Variation of Pressure distribution with respect to a and s 
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